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^ I ' Spectral analysis of the adjoint propagator in a suitable Hilbert space (and Lie algebra) 

of quantum observables in Heisenberg picture is discussed as an alternative approach to 
characterize infinite temperature dynamics of non-linear quantum many-body systems or 
quantum fields, and to provide a bridge between ergodic properties of such systems and the 
results of classical ergodic theory. We begin by reviewing some recent analytic and numerical 

V^ I results along this lines. In some cases the Heisenberg dynamics inside the subalgebra of the 

relevant quantum observables can be mapped explicitly into the (conceptually much simpler) 
Schrodinger dynamics of a single one-(or few) -dimensional quantum particle. The main body 
of the paper is concerned with an application of the proposed method in order to work out 
explicitly the general spectral measures and the time correlation functions in a quantum Ising 
spin 1/2 chain in a periodically kicking transversal magnetic field, including the results for 
the simpler autonomous case of a static magnetic field in the appropriate limit. The main 
result, being a consequence of a purely continuous non-trivial part of the spectrum, is that 
the general time-correlation functions decay to their saturation values as i~ ' . 

m 

O 

Q^^ ' §1. Introduction 

o' 

^ , Quantitative and even qualitative understanding of dynamics and ergodic prop- 

erties of interacting quantum many-body hamiltonian systems and fields (at finite (or 
infinite) temperature) is currently at its very early stage. Situation is much better, if 
p^ ! one considers bounded one or few particle systems, where ergodic theory ^' (on a clas- 

sical level) describes a variety of dynamical behaviours ranging from complete inte- 
grability being characterized by purely discrete spectrum of the appropriate evolution 
(Liouville) operator, thru ergodicity, mixing and chaos characterized by the continu- 
ed '_ ous spectrum. Mixing, which is equivalent of saying that correlation functions of an 
arbitrary pair of observables {A,B) decay in time, liint^oo{{A{t)B) — {A){B)) = 0, 
is the necessary dynamical property needed to justify the relaxation to equilibrium 
(micro-canonical) state and the laws of statistical mechanics, such as the fluctuation- 
dissipation theorem and transport laws. We know that typical classical bound few 
body system is intermediate, neither completely integrable nor fully ergodic and mix- 
ing, and hence its evolution spectrum contains both, nontrivial point spectral com- 
ponent connected to quasi-periodic dynamics in regular parts of classical phase space 
described by the KAM theorem, and the continuous spectral component connected 
to stochastic motion on chaotic components of phase space. On a quantum level. 
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bounded quantum systems of few particles have always purely discrete spectrum and 
hence their time evolution is (asymptotically) quasi-periodic, so they can never be 
truly mixing, and their dynamical properties can only in (semi) classical limit h ^ 
approach the ones of their classical counterparts. Another possibility of obtaining 
continuous evolution spectrum and truly mixing quantum behaviour is to consider 
thermodynamic limit where the size (number of degrees of freedom) of the quantum 
many-body system or field becomes infinite. It has been suggested recently^-* that 
the finite temperature dynamics of integrable quantum many-body lattices (solvable 
by quantum, inverse scattering or Bethe ansatz) is pathological (read: non-ergodic) 
from the point of view of statistical mechanics and transport phenomena, which 
may be explained by means of existence of an (infinite number of) exact conserva- 
tion laws^\ On the other hand, several numerical studies ^^'^' of high temperature 
dynamics of strongly non-integrable quantum 'many'-body systems of interacting 
particles suggested that their dynamics can indeed approach mixing behaviour in 
thermodynamic limit. However, it has been suggested in refs. ^)'^> , based on numer- 
ical results in a family of kicked fermionic lattices, that intermediate behaviour (of 
non-integrable but also non-ergodic and non-mixing quantum dynamics) may exist 
as well in thermodynamic limit in a finite range of systems' parameters. 

In ref. ^' the Heisenberg dynamics of a certain Lie subalgebra of quantum ob- 
servables equipped with a Hilbert space structure of the above mentioned kicked 
fermion model on an infinite lattice has been studied. It has been shown by numeri- 
cal algebra in operator space that the regime of the so called intermediate dynamics, 
discovered before^' by direct Schrodinger time evolution on finite lattices, exactly 
corresponds with the existence of few (in contrast to infinite for integrable case) con- 
servation laws, which are the eigenvectors (with eigenvalue 1) of the adjoint evolution 
propagator over the Hilbert space of quantum observables. 

Later on •*, a similar program has been undertaken with more analytical ap- 
proach: Dynamics of the adjoint propagator in Heiseberg picture over the two para- 
metric infinitely dimensional dynamical Lie-algebra of observables over Id quantum 
spin 1/2 chains, where the possibly time-dependent Hamiltonian can be any her- 
mitean member of the algebra, has been formally mapped onto Schrodinger dy- 
namics of a non-linear one-particle problem in 2d configuration space. Since the 
two spectral problems are shown to be equivalent, the (infinite-temperature) time 
auto-correlation functions of the spin-chains are identical to the quantum recurrence 
amplitudes of the associated one-body problem. Conceptually perhaps even more 
interesting is the result that the continuum field limit of the spin-chains corresponds 
(or maps on) to the classical limit of the associated non-linear one-particle problem 
whose dynamics can go from integrable to truly mixing and chaotic. 

In this paper we consider a related but simplified one parametric infinitely di- 
mensional dynamical Lie algebra of spin chains which has been proposed in ref. ^' and 
used to construct infinite families of conservation laws for any member of the algebra 
being interpreted as a Hamiltonian. Below we use and further develop these ideas in 
order to fully exploit the Heisenberg dynamics in the space of observables (which is 
the framework that should correspond to the classical ergodic theory ^^ around Liou- 
ville equation), and exactly compute the time-correlation functions. We consider an 
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interesting representative of this algebra, namely the Ising chain periodically kicked 
with transversal magnetic field h 

oo oo 

^Ki(i)= Y. (^^>i+i + '5r(t)K), Sr{t):= Yl Ht-mr). (M) 

where a^,p £ {x, y,z},j £ Z are the standard spin 1/2 (Pauli) operators at different 
sites j satisfying the commutation relations [o"^,o"^] = 2i6jkJ2s^prs(^j- Note that 
in the continuous-time limit r ^ (1-1) becomes an 'ordinary' Ising chain in a 
static transversal field. In sec. 2 we review some known facts ^^ about the algebra 
of Kicked Ising (KI) model. In sec. 3 we pose the spectral problem for the adjoint 
propagator in Hilbert space over the algebra ob observables and show its relation to 
correlation functions of infinite-temperature statistical mechanics. Interestingly, the 
spectral problem for the adjoint propagator can be formally interpreted in terms of a 
'quantum one-particle scattering problem' on a Id semi-infinite lattice. First, in sec. 4 
the limiting case of ordinary Ising chain in a static field is solved explicitly, and then 
in sec. 5 general results are given for the spectral measures and explicit asymptotics 
for the time- autocorrelation functions, which are shown to decay to their saturation 
values as t~^'^. However, we note a striking difference of the model in a static vs. 
periodically kicked field since the limits r ^ and t ^ oo do not commute as is 
explicitly demonstrated in case of magnetization correlation function. 

§2. Algebraic properties of Kicked Ising chain 

We start with the so-called dynaniicalLie (sub)algebra © of quantum observables 
over infinite spin chains ^' which is essentially generated by the two parts of the KI 
Hamiltonian, namely X]j CTjaj^i and ^j cr^, and which is spanned by the two infinite 
sequences of self adjoint observables C/„ and V^, n G Z, namely 

^ f ^f (^|)n-i(T|+„, n > 1, 
Un= Y { - ""h ^ = 0' (2-1) 

oo ( crJ(cr|)„_iCr|+„, n > 1, 

Vn= Y. { 1' ^ = 0' 

j=-oo ( - o-J(o-|)_„_io-J_„, n < -1, 
where ((T|)fc := Y[i=i (^j+i ^or k > 1, ((j|)o := 1, and satisfy 

[Vm, Vn] = 0, (2-2) 

One can turn the algebra 6 into the Hilbert space by defining the following (canon- 
ical) scalar product of any pair A,B £ & 

iA\B)=lun^^tTLiA^B), (2-3) 
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(tr^ is a trace over the space of finite chains of length L) with respect to which 
observables {Un,Vn',n £ Z} form an ortho-normal (ON) basis. Note that met- 
ric (2-3) is invariant with respect to the adjoint map, {adA)B = [A,B], namely 
{{ ad A'^)B\C) = {B\{ ad A)C), meaning that a hermitean observable A^ = A gener- 
ates a hermitean adjoint operator ad^ on the Hilbert space S of observables with 
respect to the metric (2-3). The canonical infinite temperature average of some ob- 
servable A £ & normalized by the 'volume' L is just (A) := {l\A). In ref. ^^ it has 
been shown that any Hamiltonian of the general form H = J2m=-m- {hmUm+ gmym) 
is integrahle since there exists analytic ©-valued function T{\) oi N = {m++m- + 1)- 
tuple of complex variables A = (Ai, . . . , Xn), |A„| < 1, commuting with the Hamil- 
tonian [H, T(A)] = for any A. From this procedure, two semi-infinite sequences 
of independent and mutually commuting conservation laws have been determined, 
the (non-trivial) charges Qk =Emi-m_ [hm{Uk+m + U^k+m) + gmiVk+m + V-k+m)], 
and the (trivial) currents Ck = Vk+i + V-k-i, for A; = 0, 1, 2 . . ., where Qo = 2H. 

Let us now turn to our KI Hamiltonian (IT), which we write as Hyii{t) = JUi — 
h6r{t)UQ, or the Floquet map factorizing into the product of kick and 'free' part 

?7k, = Texp (-i I dtH^^jitU = exp(-if ?7i) exp(if t/o) (2-4) 

where a := 2r J, and (3 := 2r/i. The key object in this paper is the adjoint propagator 
of observables in the Heisenberg picture (or the adjoint Floquet map) 

U^f = exp(-if adUo) exp(if adUi), U^fA{mT) = ulA{mT)UKi = A((m+l)r). 

(2-5) 
U^f is a unitary operator over the space of observables 6, {U^f-A\U^B) = {A\B). 
The algebra (2-2) yields a simple evaluation of the exponentials of the adjoint gen- 

erators exp(i^ adUm)A = e*2^'"Ae"*2^'", namely 

GXpyl2 adUrfijUn — C^Un ~r S^U2m—n ~r C^S.y(^Vfi—m y-in—n): 

exp(i^ adUm)Vn = C^Vn + s'?^V-n - C^Sy(Um+n-Um-n)- (2-6) 

where a shorthand notation Cj := cos 7, s^ := sin 7 is introduced. It turns out^' that 
a similar algebraic construction of T(A) and a complete set of conservation laws as 
in the Lie algebra is possible also in the corresponding Lie group. One finds 

Qk = SaCf3{Uk+i + U^k+l) - CaSp{Uk + U-k) + ^SaS/3(Vfe+l + V^k+1 " Vk-l - V-^k^l), 

(2-7) 
and the trivial currents, Ck, k >0, which are the eigenstates of the adjoint Floquet 
map with eigenvalue 1, U^fQk = Qk^ U^fCk = Ck, and generalize the known 
conservation laws for the static field ^' . Ck commute with any other element of the 
algebra & and hence span the maximal ideal 3 of the algebra &. In the following, we 
will subtract this trivial orthogonal subspace 3 and consider the Heisenberg dynamics 
on the derived algebra & = [&,&] = & — 3. 
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§3. Spectral measures and time-correlations of kicked Ising chain 

Let us now consider the spectral problem for the unitary adjoint propagator U^. 
The problem is to find, for any element A of the Hilbert space & , the corresponding 
spectral measure d^XAi^) on a unit circle ?? G [— 7r,7r), such that for a suitable test 
function f{z) one has the identity (see e.g. part III and appendix 2 of ref. ^^) 



{A\fiUi^)A) = J df^Ai^)fien. (3-1) 

The measure d^iAi'^) is composed of (a series of) delta functions for the point spectral 
component and of a continuous distribution function for the absolutely continuous 
part of the spectrum, and even of a multifractal distribution in case of a singular 
continuous spectral component. It has been shown that all the three spectral parts 
can coexist, for example for the Schrodinger problem of the quantized kicked Harper 
model ^'^^ However, we will show below that the Heisenberg propagator of KI model 
has only the trivial point spectrum t9 = corresponding to the conservation laws 
(2-7) and the absolutely continuous spectrum with the continuous spectral measure 
^a(^) ~ dfi^^' I'd) / d{} . Then, taking f{z) = z™, we write the infinite temperature 
autocorrelation function as a Fourier transformation of the spectral measure 

(A(mr)^t) = (A\[U^f]mA) = Da+ T dV^We^"'''- (3-2) 



Note that {A) = for any A £ & . Da is the time-averaged autocorrelation func- 
tion Da = limM^oo(l/2M) X]m=-Af (^('^''")^^)' o'^ ^^6 weight of the point spectral 
component and can be computed from the 'sum-rule' (putting m := in eq.(3-2)) 
Da = (^1^) — JZ-,^d'd^'^{'d). Nevertheless, Da can also be computed from the full 
set of eigenstates — orthogonalized conserved charges Q'^^ (obtained by applying the 
Gram-Schmidt orthogonalization onto the sequence Qk (2-7)), {Q'f^^Q'l) = S^i, namely 

Da = Y.\{Q',\A)\^ (3-3) 

k 

which is the essence of theorems on bounds for susceptibilities^^-*. 

Let us organize the ON-basis of & in the following way: let Eq := Uq and En be 
the triple (C/„, U^n, {Vn — V-n)/V^) for n > 1. Then general observable AG©' can 
be expanded as A = qqEq -\- Yl'^=i o-n • En ='■ a- E with one scalar and a sequence 
of vector coefficients denoted by a = (ao,ai,a2, . . .). In the basis E^ the matrix of 
adjoint map U^f can be written as a banded (3x3 block-pentadiagonal) matrix L[^, 
with a periodic structure except for small indices of rows/columns. As a consequence 
of this structure, the spectral problem for the generalized eigenfunctions ip{'d), 

U^M^) = e^^'m, (3-4) 

can be written uniquely as a 'quantum mechanical one-particle scattering problem' 
on a semi-infinite Id lattice, with the asymptotic part, for n > 2, 

A/s (Ba^n+2 - Cai^n+1 + i^aV'n + C^^n-1 + ^a V^n-2) = C^'^V'n (3-5) 
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where the 3x3 matrices Ajj , Ba ,Ca, Fa read 



Co S f. 



-V2Cf3Si3 





\/2Cf3Sf3 -\f2C(iS(i 

/O 0\ /O 0\ /c^ 

Ba = si\\ , C„ = V2c«s„ 1 , F, - ' n .2 

\0 0/ \1 0/ \0 

and with the 'scattering potential at the origin' given by the equations 

'cl-e'^ -\/2c«s« 4 o...\ 

C2« - c^e^^^ -^e^'^ -^Cf^Sf^e'' . . . U = 0. 

\/2CaSQ, \/2c^S/3e''^ -\/2c^S/3e*'' C2a - 02/36*'' -\/2cQ,Sa 0.../ 

(3-7) 
(3-7) are only equations 1,2,4 out of the first 7 rows (n < 2) of the matrix equation 
(exp(i^ adC/i) — exp(i^ adC/o)e*'')^ = which is equivalent to (3-4). The rows 3,6,7 
are included already in the asymptotic part (3-5) (for n = 1, 2 since matrices B^ and 
Ca have many zero entries), while row 5 is equivalent to row 4, and therefore they 
do not scatter the asymptotic solutions which are studied right below. 

Now we solve the asymptotic problem (3-5) with the standard ansatz leading to 

4 = A"^(A), G{\)v{\) = e'^v{\) (3-8) 

with the transfer matrix 

G{\) = ApiX^Ba - XCa + Fa + X^^Cl + X^^Bl). (3-9) 

The secular equation det(G(A) — 6X3) = reduces to a simple fourth order real 
polynomial equation in A, since det(G(A) — e^'^ls) = — ie^*''/^ sin('!9/2){a+6(A+A^^)+ 
c(A^ + A~^)}, where a = l-C2a-C2p- 3c2aC2/3 + 4 cos -;?, 6 = -2s2aS2/3, c = -2s^s^, 
with a simple solution which we write in terms of two quadratic equations 

A + A-^ = -2^^^^^±^^!^^. (3-10) 

SaSp 

Note, however that secular determinant is identically zero if the spectral parameter 
vanishes -i? = 0. Then, any function of the form (3-8) with |A| < 1 is a candidate 
for an eigenvector of U^ (or 'bound state') and it has been shown®) that linear 
combinations of three of them generally solve the boundary equations (3-7), and 
from these a complete set of (local) conservation laws (2-7) has been derived. 

In the following we thus exclude the 'trivial' eigenvalue, and fix the spectral 
parameter "i? 7^ 0. A being a solution of (3-10) means that 1/A is also a solution. We 
interpret this as conservation of the magnitude of momentum, calling A a momentum 
parameter. The corresponding eigenfunction iT(A), will be normalized for convenience, 
*' namely v{X)* • v{X) = |w(A)p = 1. Note that the eigenfunction v{X) satisfies an 

*' v{X) has been calculated explicitly using Mathematica, as all the other heavy algebraic cal- 
culations reported in this paper, and in general case is too lengthy to write down. We write it 
explicitly later in various asymptotic regimes (4-2, 5-1, 5-3). 
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interesting relation, namely v{X) ■v{X^^) = which can be proved from the following 
general property of the transfer matrix, namely G {X~^)G{X) = I3 for any A G C. 

Further, it was shown (by means of computer algebra) that the 'the scattering 
boundary condition' (3-7) can be solved for any A satisfying (3-10) with the scattering 
ansatz of an incoming wave of an amplitude 1 and outgoing wave of an amplitude 
5(A), 

V;„ = i;(A-^)A-" + 5(A)^(A)A", 



where the scattering amplitude reads 
,w • ij{X~^^ 



S{X) 



-X- 



w ■ v{X) 



with 



w 



(c; 



e-^^4'^^' 



C/3S/3) 



(3-11) 



(3-12) 




and one finds the 'unitarity condition for 
the S-matrix', |S'(A)| = 1, for either real 
or complex A. Apart from the scatter- 
ing amplitude, boundary equations (3-7) 
also determine the scalar 

iJo = V2iX-^) + SiX)v2iX). (3-13) 

From the above we learn that attenuat- 
ing wave (3-8) with real momentum pa- 
rameter A, with |A~^| < 1, cannot gener- 
ate an eigenstate of U^^ since it is always 
accompanied with exponentially growing 
wave, |A| > 1, due to the fact that 5(A) 
cannot vanish. Therefore, there is no 
point spectrum for -d 7^ 0. However, spec- 
tral parameters -d which admit complex unimodular solutions of (3-10) are in the 



Fig. 1. Spectral bands {upper - solid, and 
lower - dashed curves) for two cases: (a) 
a = 0.35,13 = 0.65, and (b) a = 1.3,/3 = 
1.1 (overlapping bands). 



absolutely continuous spectrum of the adjoint Floquet map C/^, with ^ (3- 11-3- 13) 
being the generalized eigenvectors in basis E. Writing A = e**^ in terms of quasi- 
momentum ip G [0,7r], putting it to (3-10) and solving it for i?, we obtain explicit 
forms of two continuous 'bands' ibT?(c^) with 



i9((/9) = 2arccos(cosacos/3 + sinasin/?cos(^) (mod 27r), 



(3-14) 



which go from ib2|a — (3\ for 99 = to ib2|a + (3\ for ip = tt, assuming without loss 
of generality that \a\, \f3\ < tt/2 because the full problem (3-5-3-7) is periodic in a, (3 
with period vr. Since the spectral parameter is on a unit circle the two bands overlap 
for sufficiently large kick parameters, namely if max{|a + /3|, |a — /?|} > 2^. In this 
case, one has two different quasi-momenta (p for a given fixed -d since eq.(3-10) has 
two pairs of complex unimodular solutions. See fig.l. The spectral measure of some 
observable A & & , namely fi'^{i}), is non- vanishing only inside the bands, and it may 
be rewritten in terms of quasi-momentum densities on the bands p^{ip), 



dfiAi±^iv)) = P^iv)dp, fi'Ai±Hv)) = \d^/d^\p^iip), 



(3-15) 
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where in the last equation the two terms have to be added if the bands overlap, and 
we need the 'density of states' \dip/d'd\ = \ sm{^'d) / {2sc,Sf^ sm{ip))\. 

Distributions pj^{(p) can be calculated by means of a simple-minded truncation 
of the Hilbert space 6 at n = A^ and using counting-of-states-in-a-box technique 
yielding 



PAif) 



lim 

N—>oo IT 



NKi'o + ZLi^-i^n? 



^N 



Using expressions (3-11-3-13) we obtain, writing iF := v{e^'^), S^ := S{e^^) 



(3-16) 



Pa('P) 



1 

2^ 



an[Vo 



+ S^i 



oo 



-V^-inip _j_ gtp^^imp^ 



with -& := ±'&iip). 



+ 2^ < • (^ 

n=l 

(3-17) 
We notice that p\{ip) = Pj^{^) iS A = A^ . We can finally transform the quasi- 
momentum densities back to the spectral measures, or we write the spectral decom- 
position (3-1) directly in terms of quasi-momentum integrals 



{A\f{U^f)A) 



DAf{l) + rdv9{/>+M/(e^'^('^)) +/9:iM/(e-^(^))} . (3-18) 



For example, we compute the total spectral weight of a point spectrum 
Da = {A\A) - r d^ {p\{ip) + p^iy,)} . 




(3-19) 



Formulae (3-ll-3-13,3-15,3-17,3-18) are 
very useful exact results which we can 
use with some elementary numerics in 
order to compute the spectral measures 
p'j^i'd) and time correlation functions 
{A\A(mT)). In fig. 2 we show results of 
calculation of the dynamical susceptibil- 
ity Dm, namely the time-averaged au- 
tocorrelation function of the magnetiza- 
tion M = J2j ^1 = ~Uo = m- E, m = 
(1, 0, . . .) as the function of the relative 
field strength (3/ a = h/J (for several val- 
ues of the kicking period r) . Note an in- 
teresting singularity at /i/ J = 1, or more 
generally at a = /3 (mod vr), which will 
be commented on later in sec. 5. In fig. 3 we show quasi-momentum densities 
p\{^) and spectral measures p'^{'&) for two observables, namely the magnetization 
M = —Uq and the XX-chain Hamiltonian, X = ^j(o"|(t|^-^ + a^a^j^^) = Ui + U-i, 
X = (0, 1, 1, 0, . . .), and for two different sets of parameters a, (3 (same as in fig.l, 
with and without band overlap). Moreover, one can obtain really explicit analytic 
results in two cases: (i) in the continuous-time limit r ^ of a static transversal 
magnetic field, and (ii) for asymptotically large times t = mr ^ 1. 



Fig. 



2. Dynamical susceptibility Dm vs. rel- 
ative magnetic field h/J. Different curves, 
from solid to more and more open dashed 
curve, refer to different values of r = 
0,0.25,0.5,0.75,1 (in this order). 
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Fig. 3. Quasi-momentum densities (a,b) and spectral measures (c,d) for the magnetization (solid 
curves) and XX-energy (dashed). (a,c) and (b,d) are for cases (a) and (b) of fig.l, respectively. 

§4. The limit of Ising chain in a static transversal field 



In the continuous time limit r ^ 0, t = mr, we set the energy scale by putting 
J := 1 (a = 2t, /3 = 2/it), so we are left with a single parameter h. U^f = 
1 + ir adHi + O (r^) is now an infinitesimal adjoint propagator generated with the 
Hamiltonian Hj = Ui — hUo. One can again formulate the spectral problem for the 
hermitean operator adi^i in & as (now somewhat simpler) one-particle Id scattering 
problem. Or, one merely expands general results (3-11-3-17) of the previous section 
for small r and takes the limit r ^ 0, writing the new spectral parameter of the 
hermitean operator adHi as e = -d/r. Thus one finds the bands zize{(p), 



€{ip) = 4\/l + /l2 - 2/ICOSV9 

which extend from ib4|/i — 1| to ib4|/i + 1|, the 'scattering data' 



if^= (2i{h-e''^)/€{ip), 2i{e~''^-h)/€{ip), 1 



S'^ = -l, 



(4-1) 



(4-2) 



and the density of states dyj/tie = e/(16/isinc^). Plugging all that into (3-17,3-15,3T8) 
we get explicit formulae for any observable A. Results are particularly simple (and 
perhaps physically interesting) for the magnetization M where we find 



piiiv) 



sm'^{ip) 



27r(l + /i2 -2hcosip)' 



'uie) = j^^Jm' - {h-mnh+iy - iir). 



f^M( 



(4-3) 



It is very important to note the square-root singularities of the spectral measure 
IJ,'j^{4\h — 1| + e) oc e^'^, /i'jy^(4|/i -|- 1| — e) oc e^'^, as e \ 0, since these will dominate 
the long-time behaviour of the infinite-temperature time-correlation function 



{M{t)M) =Dm + 



47r/i2 



4|?i+l| 



i\h-l\ 



decos(et)^((|)2 - {h-lY){{h+lY - (1)2). 



(4-4) 



We use an elementary asymptotics (which can be proved by complex rotation) 



iboo 



dx^J\x — a\f{x)e 



ixt 



±^/(a)exp(mt± ^isgnt)|t|-t +0 [t'^ ) , (4-5) 
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where f{x) is some analytic function, applied to both ends of the spectrum in order 
to estimate the integral (4-4) giving (again, up to O (i~^'^)) 

, ,. . |/i+l|5sin(4|/i+l|t-5)- |/i-l|5sin(4|/i-l|t+5) ^ ^ 

Further, the infinite-temperature dynamical susceptibility Dm is computed explicitly 
by integrating (3-19) the distribution (4-3) 

Z)M(/i) = l-^(max{l,|/i|})-2, (4-7) 

having the singularity at h = h^. := 1. Note that the time average Dm of {M{t)M) 
actually agrees with the r ^ limit of the general case (see fig. 2, solid curve), 
although the limits r ^ and t = mr ^ oo do not generally commute as we show 
explicitly later for the correlation function {M{t)M) itself. 

§5. Asymptotic results in general case 

Now we will establish that also spectral measures for the general case (arbitrary 
kicking parameters a,/3) have square-root singularities (see fig.3) and thus lead to 
^-3/2 (jggay of correlations. We will first only consider contributions from the upper 
spectral band +'d{(p)^ while the contribution of the lower band is obtained simply 
by replacing A by A'^ (or a by a*). Below we assume that a > 0,/3 > 0, so |a — 
P\ ^ |a + /3|, while other cases can be obtained with trivial modifications. Let 
us first expand the band around the minimum for small quasi-momentum ip = e, 
i?(e) = 2|a-/3| + i(saS/3/s|a-/3|)e^ + (e^). Then the scattering data (3-8,3-ll-3-13) 
are expanded explicitly in leading two orders 

„2isgn(a-/3)a _ i 

i^ = i-h h «gn(a-/3)^) - Is ^,,|,_,|_^ (1, 1, 0) + O {e') , (5-1) 

When we evaluate the quasi-momentum densities (3-17) we find that expression 
inside |.| vanishes to order O (1) so we have p\{6) = Ka£^ + O (e^) with coefficient 

Ka := — 

27r 



ao + Va„-|^^(l,l,0)+ne^^g°(°-'^)'^(-l,l,isgn(a-/3)^/2)| 

(5-2) 



Sa-I3 



Similarly we find, expanding around the other end (the maximum) of the band, 

V? = vr — e, ■d{iT — e) = 2\a + /3| — ^^{^aSp/ s\a+p\)^'^ + C* (e^), the scattering data 



1 e^^" - 1 
v" — ' " " " ' 



and the quasi-momentum densities p\{t^ — e) = Las"^ + O (e^) with 



... - 2 

La ■= 7— ao- 

ivr Sr 



00 ^ ^ 

I^ + V(-l)"a„ • ^(1, 1, 0) + ne-^(-l, 1, -zV2) 



(5-4) 



Kicked Ising chain 
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Transforming to spectral variable -d and multiplying by the density of states (3-15), 
which has a simple form oc j"!? — 2\a it /3|p^ around both respective ends of the 
spectral band, we obtain explicit square-root singularities of the spectral measure 



3l + «) = i 



fi'A2\a + /3\ - 



Sa-I3 



SaSp 

Sa+13 



SaSf3 



KaV^ + O 



and similarly for the lower band by replacing A by A' . Since the spectral measure 
n'j^i'd) is a smooth- function on a complex unit-circle, except for four singularities at 
the four band edges ib2a it 2/3, the asymptotic approximation to the integral (3-2) is 
dominated by the four terms which are computed using asymptotics (4-5) (see fig. 4) 



Sa-P 



{A{mr)A^)^DA + ^{ 



+ 



SaSp 
Sa+P 



SaSf3 






with Tj := (37r/4)sgnm. 



m\ 2 
(5-5) 



In case of a hermitean operator A = A' the formula simplifies and then, of course, 
time-correlations are real and symmetric {A{mT)A) = {A{—mT)A) = {A(mT)A)* . 
For example, the asymptotics for magnetization is (for m > 0, again upto O (m~^''^)) 



{M{mT)M) ^Dm + 



1 



■KS 



sin(2|a-F/3|m-|) s\n{2\a-(3\m+\) 



5a+/3 



1/2 



11/2 



m 



\Sa-fi\ 

(5-6) 
which does not converge to (4-6) if we let r = a = (3/ J -^ while keeping t = mr 
large meaning very explicitly that the limits r ^ and t ^ oo do not commute. 

Note that all the quantities computed above are non-smooth functions of pa- 
rameters a,f3 on the line a = P (mod n) (e.g, see fig. 2) since the band minimum 
^min = 2\a — P\ is non-smooth. For a = (3 the two bands touch (at the point t? = 
where /^^(t?) then becomes a smooth function), and the square-root singularities at 
±2\a — P\ = dissapear, and so should also the two terms with Ka, Kj^t in (5-5). 

§6. Conclusion 

The problem of infinite-temperature time-correlation functions in a family of 
Ising spin 1/2 chains periodically kicked with transversal field, which has been for- 
mulated in terms of a spectral problem for the adjoint propagator over a certain 
subspace of observables in Heisenberg picture, has been solved using methods (and 
terminology) of a single-particle quantum scattering on a semi-infinite tight-binding 
lattice. It has been shown that time-autocorrelation function generally decays as 
t~^" to its saturation value which is, due to integrability, generally different from 
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0.02 



100 




100 



Fig. 4. The time-correlation functions (average subtracted) of (a) magnetization M for a — 0.35, 
l3 = 0.65, and of (b) XX-energy X for a = 1.3,/3 = 1.1. Dots are exact (numerical) results while 
thin solid curves are asymptotic formulae (5-5, 5-6). 



the squared canonical average. This may be interpreted in terms of a relaxation to a 
non-unique equihbrium statistical steady-state. Furthermore, it excludes the possi- 
bility of (quasi) periodic motion, which is a qualitatively different situation than the 
one we encounter in few-body classical (or quantum) integrable systems. Note that 
such behaviour is drastically different from dynamics of some 'trivially' integrable 
quantum many-body lattices, such as XX or Ising chain without external magnetic 
field (e.g., put /3 = /i = in the results above) where the continuous spectrum of 
adjoint dynamics collapses to a point and one recovers periodic time-correlations. It 
is an open challenge whether our approach can be extended to more 'sophisticated' 
integrable quantum lattices, such as the general XYZ-chain or the Hubbard model, 
perhaps within a formalism of quantum inverse scattering. 

I hope that these results may find some interesting application, and may help 
to stimulate some development of ergodic theory of quantum many-body systems. 
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